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Abstract. In this article we investigate the existence of a solution to a semi-linear, 
elliptic, partial differential equation with distributional coefficients and data. The prob- 
lem we consider is a generalization of the Lichnerowicz equation that one encounters in 
studying the constraint equations in general relativity. Our method for solving this prob- 
lem consists of solving a net of regularized, semi-linear problems with data obtained by 
smoothing the original, distributional coefficients. In order to solve these regularized 
problems, we develop a priori i°°-bounds and sub- and super-solutions to apply a fixed 
point argument. We then show that the net of solutions obtained through this process sat- 
isfies certain decay estimates by determining estimates for the sub- and super-solutions 
and utilizing classical, a priori elliptic estimates. The estimates for this net of solutions 
allow us to regard this collection of functions as a solution in a Colombeau-type algebra. 
We motivate this Colombeau algebra framework by first solving an ill-posed critical ex- 
ponent problem. To solve this ill -posed problem, we use a collection of smooth, "approx- 
imating" problems and then use the resulting sequence of solutions and a compactness 
argument to obtain a solution to the original problem. This approach is modeled after 
the more general Colombeau framework that we develop, and it conveys the potential 
that solutions in these abstract spaces have for obtaining classical solutions to ill-posed 
non-linear problems with irregular data. 
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1. Introduction: Semilinear Problems and Critical Nonlinearities 

In this paper we consider a family of elliptic, semilinear Dirichlet problems that are of 
the form 

N K 

- Di(a ij D jU ) + ^Vu n * = in Q, (1.1) 

i,j=l i=l 

u\an = P, (1.2) 

where a l \ b % and p are potentially distributional and n, E Z for each i These problems 
are a generalization of the the Lichnerowicz Equation that appears in the study of the 
constraint equations of General Relativity. The need to understand an equation of this 
form with rough data arises if one attempts to study equations such as the Lichnerowicz 
equation when the metric of the embedded hypersurface is not smooth. From a physical 
perspective, such problems are interesting because distributional metrics correspond to 
the initial data for physically plausible spacetimes generated by strings and gravitational 
waves (H [6]|. From a mathematical point of view, developing results for solutions to 
the Lichnerowicz equation under such low regularity conditions is of interest in that it 
extends the current "rough metric" existence theory, as described in |[8l [TQl |H |2l . The 
hope is to eventually extend the solution theory to cover rough data examples such as 
those studied by Maxwell in [fTTTl . 

The main contributions of this article are an existence result for (11.11) in a Colombeau- 
type algebra, and an existence result in W 1,2 (Vt) for an ill-posed, critical exponent prob- 
lem of the form 

- Am + au 5 + bu { = in fi, (1.3) 

u\an = p, 

where 1 < i < 4 is in N, Q C R 3 , a E L p (Vt) and b E L q with § < p < q. The 
framework we use to prove existence for (11.11) consists of embedding the singular data 
and coefficients into a Colombeau-type algebra so that multiplication of the distributional 
coefficients is well-defined. To solve (11.31) . we do not explicitly require the Colombeau 
machinery that we develop to solve (|l.ll) . but we use similar ideas to produce a sequence 
of functions that converge to a solution of (11.31 ) in W l,2 (VL). 

The Colombeau solution framework for this paper is based mainly on the ideas found 
in lfi~2l . Here we extend the work done by Mitrovic and Pilipovic in lfT2ll to include 
a certain collection of semilinear problems. While Pilipovic and Scarpalezos solved a 
divergent type, quasilinear problem in a Colombeau type algebra in lfT3l . the class of 
nonlinear problems we consider here does not fit naturally into that framework. Here we 
provide a solution method that is distinct from those posed in |[T3l and [fl2ll that is better 
suited for the class of semilinear problems that we are interested in solving. The set up 
of our problem is completely similar to the set-up in lfi~2ll : given the semilinear Dirichlet 
problem in (11.11 ), we consider the family of problems 

P e (x,D)u e = f e {x,u e ) onfi, (1.4) 

u e \an = Pe, 

where f e , h e , and P e (x, D) are obtained by convolving the data and coefficients of (11.11) 
with a certain mollifier. Thus a solution to the problem in a Colombeau algebra is a net 
of solutions to the above family satisfying certain decay estimates in e. This is discussed 
in detail in Sections [3T21 and [3741 This basic concept underlies both the solution process 
in our paper and in |[T2ll and [fT3l . However it is our solution process in the Colombeau 
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algebra that is quite distinct from that laid out in lfT2l . where the authors used linear 
elliptic theory to determine a family of solutions and then classical elliptic, a priori 
estimates to prove certain decay estimates. Most notably, the authors developed a precise 
maximum principle-type argument necessary to obtain the decay estimates required to 
find a solution. Our strategy for solving (11.11) differs in a number of ways. First, in 
Section I5TT1 we develop a family of a priori L°° bounds to the family of problems (11.41) . 
Then in Section [521 we show that these estimates determine sub- and super- solutions to 
(11.41) . We then employ the method of sub- and super-solutions in Section l4Tl to determine 
a family of solutions. Finally, e-decay estimates on the sub- and super-solutions are 
established in Section [521 an d in Section [6] these estimates are used in conjunction with 
the a priori estimates in Section 13.11 to prove the necessary e-decay estimates on our 
family of solutions. 

This paper can be broken down into two distinct, but related parts. The first part is 
dedicated to solving (|1.3I) . Our solution to this problem does not explicitly require the 
techniques that we develop to solve problems with distributional data in Colombeau al- 
gebras and only relies on standard elliptic PDE theory. However, the ideas that we use 
to solve the problem are closely related: we obtain a solution by solving a family of 
problems similar to (11.41) and then show that these solutions converge to a function in 
W 1,2 (Vt). Therefore, we present our existence result for (11.31) first to convey the benefit 
that the more general Colombeau solution strategy has, not only for solving problems in 
the Colombeau Algebra, but also for obtaining solutions in more classical spaces. The 
remainder of the paper is dedicated to developing the Colombeau framework described 
in the preceding paragraph. This consists of defining an algebra appropriate for a Dirich- 
let problem and properly defining a semilinear elliptic problem in the algebra. Once a 
well-posed elliptic problem in the Colombeau algebra has been formed, we discuss the 
conditions under which the problem has a solution in the algebra and finally, describe 
how to translate a given problem of the form (11.11) into a problem that can be solved in 
the algebra. It should be noted that while the intention is to find solutions to (11.11 ), the 
main result pertaining to Colombeau algebras in this paper is Theorem 14 .11 which is the 
main solution result for semilinear problems in our particular Colombeau algebra. 

Outline of the paper. The remainder of the paper is structured as follows: In Section [2] 
we motivate this article by proving the existence of a solution to (11.31) . In Section [3] 
we state a number of preliminary results and develop the technical tools required to 
solve (ll.l|) . Among these tools and results are the explicit a priori estimates found in |[T2ll 
and a description of the Colombeau framework in which the coefficients and data will 
be embedded. Then in Section [4] we state the main existence result in Theorem 14. 11 give 
a statement and proof of the method of sub- and super solutions in Theorem 14.31 and 
then give an outline of the method of proof of Theorem 14.11 Following our discussion 
of elliptic problems in Colombeau algebras, we discuss a method to embed (11.11) into 
the algebra to apply our Colombeau existence theory. The remainder of the paper is 
dedicated to developing the tools to prove Theorem 14.11 In Section [5] we determine a 
priori L°° bounds of solutions to our semilinear problem and a net of sub- and super- 
solutions satisfying explicit e-decay estimates. Finally, in Section[6]we utilize the results 
from Section \5\ to prove the main result outlined in Section [4] 

2. SOLUTION CONSTRUCTION USING A SEQUENCE OF APPROXIMATE PROBLEMS 

If O C M 3 , the Sobolev embedding theorem tells us that W l,2 {Vt) will compactly 
embed into L p (Vt) for 1 < p < 6 and continuously embed for 1 < p < 6. Given functions 
u,v E W l,2 (Vt), this upper bound on p places a constraint on the values of i that allow for 
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the product u l v to be integrable. In particular, Sobolev embedding and standard Holder 
inequalities imply that this product will be integrable for arbitrary elements of W l,2 (Vl) 
only if 1 < i < 5. More generally, if a G the term au 5 v will also be integrable. 

However, if a is an unbounded function in L p {Vt) for some p > 1, then this product is not 
necessarily integrable without some sort of a priori bounds on a, u, and v. Therefore, the 
following problem does not have a well-defined weak formulation in W 1,2 (£l): 

-Au + cm 5 + bu { = inQ, (2.1) 

u\an = P, 

where 1 < i < 4 is in N, p G W 1 ' 2 ^'), a G L P (VL'), b G L q (VL') for § < p < q and 
where CC O' are bounded domains in M 3 . 

The objective of this section is to find a solution to the above problem. In order to solve 
(12.11 ), we solve a sequence of approximate, smooth problems and use a compactness 
argument to obtain a convergent subsequence. We first define necessary notation and 
then present the statements of two theorems that will be necessary for our discussion in 
this section. Then we prove the existence of a solution to (|2.1I) . Finally, we show that if a 
solution exists, then under certain conditions we can construct a net of problems whose 
solutions converge to the given solution. 

2.1. Overview of Spaces and Results for the Critical Exponent Problem. For the 

remainder of the paper, for a fixed domain Q C R n , we denote the standard Sobolev 
norms on Q by 



H\lp = (J \u(x)\ p dx ) , (2.2) 
( k 

\u\\ W k, v = y \\d i u\ 



I LP 



Furthermore, let 

ess sup u = u, (2.3) 
ess inf u = u. 

In our subsequent work we will also require regularity conditions on the domain Vt and 
its boundary. Therefore, we will need the following definition taken from Q : 

Definition 2.1. A bounded domain Q C W 1 and its boundary are of class C k,a , < a < 
1, if for each x G dVL there is a ball B(xq) and a one-to-one mapping ^ of B onto 
D CW 1 such that: 

(1) $(5nfi) C M™, 

(2Jf(Bn dQ) C dRl, 

(3) ^ G C k ' a (B), G C k ' a (D). 

We say that a domain is of class C°° if for a fixed < a < 1 it is of class C k ' a 
for each k G N. Additionally, for this section and the next we will require the following 
Theorem and Proposition: 

Theorem 2.2. Suppose Q C E n is a C°° domain and assume f : Q x M + — > M. is in 

C°°(Q x M + ) and p G C°°(fl). Let L be an elliptic operator of the form 

Lu = —Di(a ij Dju) + cu, and a ij , c G C°°(n). (2.4) 
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Suppose that there exist sub- and super-solutions u_ : Q — > R and u + : Q — > R such 
that the following hold: 

(1) u-,u+eC°°(n), __ 

(2) < u-(x) < u+(x) Wx G Q. 

Then there exists a solution u G C 00 (r2) to 

Lu — f(x,u) onQ, (2.5) 

u\ d n = P, (2.6) 

such thatu-(x) < u(x) < u + (x). 

Proposition 2.3. Let ube a solution to a semilinear equation of the form 

N K 

- Di(a ij D jU ) + hiyni = in ( 2 - 7 ) 

i,j i=l 

u \an — Pi p( x ) > on dfl 

where a l \b l and p G C°°(Q). Suppose that the semilinear operator in (|2.7I) has the 
property that rii > for some 1 < i < K. Let uk be the largest positive exponent and 
suppose that bx (x) > in fl. Additionally, assume that one of the following two cases 
holds: 

(1) fii < for some 1 <i < K andifrii = min{nj : rii < 0}, (2.8) 

then b\ (x) < in Q. 

(2) uk is odd and < rii for alll < i < K. (2.9) 

If case (12.81 ) holds, define 

a' = sup \ sup bi(x)b ni <0 V6e(0,c)l, (2.10) 
and let a = min{a', inf xg gn p(x)}. If case (|2.8I) or case (12.91 ) holds, define 

(3'= inf (V inf bi{x)b n > >0 V6g(c,oo)1, (2.11) 

cGK I ^— - ' x£Q 

and let (3 = m&x{(3', sup x( z dn p(x)}. 

Under these assumptions and definitions, if case (12.81) holds, then < a < u < (3. 
Otherwise, if case (12.91 ) holds, then —f3'<u<(3. 

For a more detailed statement of Theorem 12 .21 and its proof, see Section |4TT] The proof 
of Proposition [23] can be found in Section [51 Now that we have all of the tools we need, 
we shall now prove the existence of a solution to a problem of the form (12.11) . 

2.2. Existence of a Solution to an Ill-Posed Critical Exponent Problem. For the fol- 
lowing discussion, let Vt C R 3 be a closed and bounded domain of C7°°-class. We assume 
that H CC O' C R 3 , with Vl' open and bounded. Here we prove the existence of a solu- 
tion to the problem 

-Am + au 5 + bu { = in fi, (2.12) 

u\an = p, 
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where 1 < i < 4 is in N, 

aeL p {Q'), b G L q (Q') fl L°°(Q'), ^ < p < q < oo, peW 1 ' 2 ^'), (2.13) 

and 

a>0, 6<0, and p > 0. (2.14) 

Proposition 2.4. The semilinear problem (12.121) /las a solution u G jy 1,2 (fi) z/a, b, and p 
satisfy the conditions in (12.131 ) and (12.141 ). 

Proof. To determine a solution to (12.121) . we consider the sequence of solutions to the 
approximate problems 

-Au n + a n (u n f + b n (u n y = infi, (2.15) 

u n \an = Pn, 

where a n = a*<p n , b n = b*<fi n , and p n = p*<fi n and n = n 3 0(nx) is a positive mollifier 
where J dx = 1. Given that is a positive mollifier, it is clear that for each n G N, 

a„ > 0, 6„ < and p n > 0. 

We first verify that the sequence of problems (12.151) has a solution for each n. To do this, 
we will utilize Theorem |2.2| and Proposition [231 Let f3' n and f3 n have the same properties 
as f5' and (3 in Proposition 12.31 for the sequence of problems (12. 1 5b . Then using the 
notation in Proposition |2.3l we can write explicit expressions for f3 for (12.121 ), f3' n and f3 n . 
It is not hard to show that 

= max 
and 

Pn = \ J and 0n = max {ft p n } . 

V a nJ 

By Proposition 12.31 f3' n and f3 n are a priori bounds for the approximate problems. Fur- 
thermore, it is not difficult to see that for each n G N that — (3' n and f3 n are sub- and 
super- solutions for ( 12.151) . See Section [531 and Theorem 14.31 for more details. Therefore, 
given that p n , a n , b n G C°°(fi) for each n, we have that u n G C°°(£l) is a solution to 
(12.151 ) for each n by Theorem |2.21 

Now observe that for each n G N, (3 n < (3, which follows from the fact that 

-b n {x) = [ (-%))0 n (x - y) dy < [ (-b)4> n (x -y) = -b, (2.16) 




and a n (x) > a, which is verified by a similar calculation. Therefore, by standard L p 
elliptic regularity theory 

||«n||w2>p <C(|| - a n (u n ) 5 - b n (u n y\\ LP + ||-u n || LP ) (2.17) 

<C(f3l\\a n \\ LP + f3 l n \\b n \\ LP + f3 n ) < M < oo, 

where M is independent of n given that f3 n < (3, a n — > a in L p , b n — > b in L q and 
p < q. Because p > | and f2 is of C°°-class, W 2,P (VL) embeds compactly into W 1,2 (Vt). 
Therefore, there exists a convergent subsequence u nj — > u in W l,2 (Q). We claim now 
that u satisfies the following two properties: 

(1) —ft < u < (3 almost everywhere, 

(2) u solves (l2~12l) . 
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The inequality —f3 < u < (3 a.e. follows from the fact the u n . — > u in W 1,2 (Q) and 

~P < -Pn 3 < -P' n . < U nj < n . < (3. 

Indeed, if we assume that u > (3 on some set of nonzero measure, then for some n the 
set A n = {x E Vl : u(x) > (3 + ^} has positive measure. Then for all j E N, we have 
that 

/ \u ni - u\ 2 dx> j \u n . - u\ 2 dx > —fj,(A n ) > 0. 

But this clearly contradicts the fact that u Uj — > u in W /1,2 (fi). A similar argument shows 
that u > —/3, a.e in fi. 

Finally, we want to show that u solves (|2.12l) . Let e > 0. Then we have that 

J (Vu-Vv + au b v + bu i v) dx = J (Vit -Vv + au 5 v + bu i v) dx (2.18) 



j (Wu nj ■ Vv + a nj {u nj fv + b nj (u nj ) l v) dx 



given that u nj solves (|2.15l) . Then expanding the second line of the above equation we 
find that 



Vu ■ Vv + au v + bu l v dx 



(2.19) 



< / | Vn ■ Vf — Vu nj ■ Vv | dx + \au v — a n . (u n .) v\ dx 



- / \bu l v - b n .(u n .) % v\ dx 



(2.20) 



< / I Vu • Vv — Vu n . ■ Vv\ dx + \au 5 v — a(u n .) 5 v I dx 



+ J \a(u nj ) 5 v - a n .(u nj ) 5 v\ dx + / \bu l v - b{u nj ) l v\ dx 
+ / \b(u n yv - b n Au nj ) l v\ dx. 



(2.21) 



Every term in (12.211) tends to given that u Tl] — > u in W X ' 2 {Q), a nj — > a in L p (fl), 
b n . — > b in L q (Q) and — f3 < u < (3. To show that the expression 



\au 5 v — a(u n .) 5 v\ dx — > 0, 



we can use a power series expansion to obtain 



\au 5 v — a(u n .) 5 v I dx 



i=l ^ ' 

5 f 

<C(5,/3 4 )^^ / | (m — u nj )av\ dx 
i=i J 

< C\\a\\Lp\\u — u n .\\ 2 P \\v\\ 2 P 



dx 



< C\\a\\ L p\\u — u ni \\ w i,2\\v\ 



w 1 
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where the last inequality follows from the fact that W ' (ft) embeds into Lp- 1 (fi) if 
p > |. Finally, by definition of the trace Tu, we have that, 

Tu = lim u n .\ d Q = lim p n? = p, 

where the limit is taken in L 2 (dVl). Therefore u — p G Wq' 2 (£1) and so u solves (12. 1 2b . 

□ 

2.3. Convergence of Approximate Solutions to an Existing Solution. In this section 
we again assume Q C M 3 is of C°°(^)-class and that Q CC Q', with Q' C R 3 open and 
bounded. We also consider the same semilinear problem as in the previous section: 

-Am + au 5 + bu i = in Q, (2.22) 

u\an = P, 

and assume that a G L P (Q'), b G L q (Q!), § < p, q < oo, p G W 1 ' 2 ^') and 

a > 0, b < 0, and p > 0. 

Now we assume that a solution u G Jy 1,2 (f2) to (12.221) exists, and we consider the con- 
vergence of solutions to the following net of approximate problems to u: 

-Au e + a e {u e f + b e (u € y = infi, 

u e \dn = Pe, (2.23) 

where a e = a * <fi e , b e = b * <fi e , p e = p * <p e , 4> e (x) = \4>{-) and is a positive mollifier 
such that J cj)(x) dx = 1. Again, given that is positive, for each e < 1 

d € > 0, b e < 0, and p e > 0. 
Because u is a solution to (|2.22l) . Proposition 12 . 3 1 implies that 

-0<-(-|) <M<max|p, (^-1^ }>=/3. 
Similarly, for each e G (0, 1) we have that 

i 



-/?e < M e < max |p e , l^-^J J = 

Given that e is a positive mollifier, (3 e < f3 for all e G (0, 1). This follows from (12.161) . 
Therefore -(3 <u e < (3 for all e G (0, 1). 

Proposition 2.5. If a solution u to (12.221) exists in W 1,2 (Q), then the approximate solu- 
tions u e of (12.231) converge to u in W 1,2 (fl) provided that 

maxjp, (^-0 J«l. (2.24) 

Proof. Define u e to be the solution to the following auxiliary problem: 

-Au e = -a e u 5 - b e u { in Q, (2.25) 
u £ \an = Pe- 

Note that u e exists by standard linear, elliptic solution theory because of the bounds on u 
and the assumption that p, q > | implies that 

-a e u 5 -bj G H-\Q). 
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Now consider the following string of inequalities: 

\\u — UtWw 1 ' 2 — ||w — M e ||v^i,2 + \\u e — UeWw 1 - 2 (2.26) 
<C {\\(a e -a)u b + (b e -b)u% H ~, 

+C\\a e (u e ) 5 + b e (u e y - a e u b - b e u A \\ H -i) . 

We observe that 

||(a e - a)u 5 + (b e - byWn-i < (3 5 \\a e - a\\ LP + (5 l \\b t - b\\ L * ->• 0, (2.27) 

for any 1 < p, q < oo. 

Furthermore, we may rewrite the second term in the inequality using a power series 
expansion to obtain the following string of inequalities: 

\\a t (u e f + b t (u e y - a € u 5 - b e u l \\ H -i (2.28) 

< \\a e J2 (^j(«> e ) 5 - j (u-u e y\\ H -i + \\b e j2^(u e y-^u-u e y\\ H -,. 

Given that \u — u e \ < 2/3 a.e. , for 1 < p' < p, we have 

5 /^x 5 



i=i j=i 

<2 4 /3 4 C(5,fi)||a e (u-u e )|| Lj/ 



LP 



< C7(^, 5)/9 ||a e || iP ||u — w e || pp > 

Lp-p' 

4| 



<C(Q,5)/3 ||a e || LP (o)||M - mJvk 1 . 2 , (2.29) 
provided that p > |. Similarly, if g > |, we have that 

H 6e £ (^M~ j ( u - u *) S \\H-i < C(4,n) / 9 i - 1 ||6 e || w ||«-tt e || W ri. a . (2.30) 

Therefore, equations (12.261) - (12.301) imply that 

||w — Mell^ 1 - 2 < /3 4 \\a e — cl\\lp + /3 4 \\h — 6||l<? + C(f2, S^HajLpHw — WeHtf/ 1 ' 2 

+ C7(4, nj^-^feelli, ||U - Uell^i.a. (2.31) 

Given that \\a e — a\\zp — > and ||6 e — b\\ L g — )■ 0, if 

/3 4 ||a|| LP « 1, and /^l&lb « 1, 
we will have that \\u — u e \\ w i,2 — > 0. But (12.241) implies the above condition provided 



that 



max ^ /,.[--] } <:<: I. 



is sufficiently small. □ 

Now that we have demonstrated the potential that sequences of solutions to approxi- 
mate problems have in solving classical semilinear equations, we begin to develop our 
Colombeau Algebra framework. Within this framework, one embeds a problem with 
rough data into the algebra and obtains a net of solutions. If the net satisfies certain de- 
cay estimates, it is regarded as a member of the algebra. Then ideally one can relate the 
net of solutions to a classical solution using techniques like those discussed in section |2T2l 
andE: 
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3. Preliminary Material: Holder Spaces and Colombeau Algebras 

We now begin to develop the Colombeau Algebra framework that will be used to solve 
( ll.lt . We first define Holder Spaces and state precise versions of the classical Schauder 
estimates given in Ifl2ll . The definition of the Colombeau Algebra in which we will 
be working and these classical elliptic regularity estimates make these spaces the most 
natural choice in which to do our analysis. Therefore we will work almost exclusively 
with Holder spaces for the remainder of the paper. Following our discussion of function 
spaces, we define the Colombeau algebra in which we will work and then formulate an 
elliptic, semilinear problem in this space. 

3.1. Function Spaces and Norms. In this paper we will make frequent use of Schauder 
estimates on Holder spaces defined on an open set f2 C W 1 . Here we give notation for 
the Holder norms and then state the regularity estimates that will be used. 

All notation and results are taken from (5]|. Assume that f2 C W 1 is open, connected 
and bounded. Then define the following norms and seminorms: 

\u{x) -u{y)\ 

[u\a;Q = SUp - , (3.1) 

x,yGQ \X y\ 

x^y 

Mfc,o ; n = sup sup \D p u\, (3.2) 
\p\=k xen 

[u]k,a-,n = sup [D p u)a;a, (3.3) 
\p\=k 



k 



\\ u \\c*(n) = \ u \k-,n = 2_j[u]j,o;n, (3.4) 

j=0 

\\ U \\c k ' a (n) = \ u \k,a;Q = \u\k;Sl + [w]fc, a; fi. (3.5) 

We interpret C k ' a (Vl) as the subspace of functions / G C k (Vl) such that is a-Holder 
continuous. Also, we view the subspace C k,a (VL) as the subspace of functions / G C k (VL) 
such that p k > is locally a— Holder continuous (over compact sets K CC Vt). 
Now we consider the equation 

Lu = a tj DijU + UuDiU + cu = f in VI, (3.6) 

u = p on dil, (3.7) 

where L is a strictly elliptic operator satisfying 

a ij = a ji and a ij (x)Uj > Mt\ 2 , ^ e Q, ^1". 

The following regularity theorems can be found in [5] and [12]. See JH for proofs. Note 
that the constant C in the following theorems has no dependence on A or A. 

Theorem 3.1. Assume that Q is a C 2,a -class domain in W 1 and that u G C 2,a (Q) is a 
solution (|3.6t . where f G C a (Q) and p G C 2,a (Q). Additionally assume that 

\ a ' lJ \0,a;Q, \b l \o,a;Q, |c|o,a;Q — A. 

Then there exists C > such that 

\u\2,a;U < C (^j (\u\ ;Q + \p\^,a;Sl + |/|o,a ; n)- 
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This theorem can then be extended to higher order derivatives if one inserts the deriva- 
tive of a solution u into (13.61 ) rearranges the equation and repeatedly applies Theorem 13.11 
See [12] for details. We summarize this result in the next theorem. 

Theorem 3.2. Let Vt be a C k+2 ' a -class domain and u G C 2 Vt fl C°(Q) be a solution 
of (13.61) . where f G C k ' a (Q) and p G C k+2 > a (Q). Additionally assume that 

Thenue C k+2 ' a ' n (Tl) and 

/ A \ 3(fc+l) 

|w|fe+2,a!;n < C + [ — J {\ u W,n + \p\k+2,a;Q + \f\k,a;U), 

where C is the constant from Theorem \3.1\ 

3.2. Colombeau Algebras. Now that we have defined the basic function spaces that we 
will be working with and stated the regularity theorems that will be required to obtain 
necessary decay estimates, we are ready to define the Colombeau algebra with which we 
will be working and formulate our problem in this algebra. 

Let V be a topological vector space whose topology is given by an increasing family 
of seminorms fi k . That is, for u G V, Hi(u) < Hj(u) ifi< j. Then letting / = (0, 1], we 
define the following: 

Ey = (V) 1 where u G Ey is a net (u e ) of elements in V with e G (0, 1], (3.8) 

£ My = {( Ue ) e E v | VA; G N 3a G R : fi k (u t ) = 0(e a ) as e 0}, (3.9) 

My = {(O G E VM | Vfc G N Va G R : fi k (u € ) = 0(e a ) as e 0}. (3.10) 

Then the polynomial generalized extension of V is formed by considering the quotient 
Qy = Em,v/Nv- Let's give a few examples, following [fT2l l6l. 

Definition 3.3. IfV = C, r G C, /^(r) = |r|, one obtains C, the ring of generalized 
constants. 

Definition 3.4. Let Q C M. n be an open set, U k C C Q an exhaustive sequence of compact 
sets and a G Nq a multi-index. Then if 

V = C7°°(n), / G C°°(Q), Mfe (/) = supd^/l : x G U k , \a\ < k}, 

one obtains Q s (fl), the simplified Colombeau Algebra. 

Definition 3.5. IfV = C°°(f2), where flcK" /.v bounded and 

f i k (f)=sup{\D a f\:\a\<k, x G Q}, 

we denote the generalized extension by Q (Q). The set E M Coc «?> will be denoted by Em (^) 
and be referred to as the space of moderate elements. The setN Caa (m will be denoted by 
Af(Q) and will be referred to as the space of null elements. 

Both Q S (Q) and C were developed by Colombeau and laid the basis for the more 
general construction described in (I3.8I) - (I3.10I) . See for more details. As in (12) , for the 
purposes of this paper we are concerned with Q(Vt) given that we are interested in solving 
the Dirichlet problem and require a well-defined boundary value. If (u e ) G £m(Q) is a 
representative of an element u G G(VL), we shall write u = [(u e )] to indicate that u is 
the equivalence class of (u € ). At times we will drop the parentheses and simply write 
[u e ]. Addition and multiplication of elements in Q{Q) is defined in terms of addition and 
multiplication of representatives. That is, if u = [(u e )} and v = [{v e )], then uv = [(u e v e )] 
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and u + v = [(u e + v e )]. Derivations are defined for u = [(u e )] E Q(Vi) by d Xi u = 
[(d Xi u e )}. 

Theorem 3.6. With the above definitions of addition, multiplication and differentiation, 
is a associative, commutative, differential algebra. 

Proof. This follows from the fact component-wise addition, multiplication, and differ- 
entiation makes V 1 = (C 00 ^)) 1 into a differential algebra. By design, £m{&) is a the 
largest sub-algebra of (C 00 ^)) 1 that contains Af(Q) as an ideal. Therefore Q{Vt) is a 
differential algebra as well. See |6]|. □ 

Now that we have given the basic definition of a Colombeau algebra, we can discuss 
how distributions can be embedded into a space of this type. 

3.3. Embedding Schwartz Distributions into Colombeau Algebras. While the alge- 
bras defined above are somewhat unwieldy, these spaces are well suited for analyzing 
problems with distributional data. The primary reason for this is that the Schwartz distri- 
butions V'iVt) can be linearly embedded into them. This allows one to define an extrinsic 
notion of distributional multiplication that is consistent with the pointwise product of C°° 
functions. Here we discuss a method of embedding T>'(Q) into Q s {Vt) in the event that Q 
is a bounded open subset of MJ 1 . Recall that G s (£l) was defined in (13 .4b ■ 

The Schwartz distributions on an open set f2 C M n , denoted V'(Q), are defined to be 
the dual of V(Q), the space of C°° functions with support contained in fi. For a given 
if E V(Q) and T E P'(O), the action of T on <p will be denoted by (T, tp). We now 
define an embedding of V(Vt) into Q s {Vt) and state some of its properties without proof. 

We begin by letting ij) E SiW 1 ), the Schwartz functions, be a function such that tp = 1 
on some neighborhood of 0. Then define E SiW 1 ) by = J 7-1 [tp], the inverse Fourier 
transform of ip. It is easy to see that 

/ <f)dx = l and / x a <p dx = V|a| > 1. (3.11) 

Let (fi t = ^-0(f ) be the usual mollifier. 

The properties of specified in (13.1 II) are extremely important. In particular, if C 
W 1 is open and bounded and ^ = 1 in a neighborhood of f2, this choice of as a mollifier 
makes following map well-defined: 

% : V\Q) G s (f2), (3.12) 

*(«) = (((e«)*0 e )ifl)- 

We summarize some properties of this map in the following theorem: 

Theorem 3.7. Let Q C M n be a bounded open set and let £ = 1 in a neighborhood ofQ. 
Define the map 

cr :C°°(f2) — y Q s (Jtl), (3.13) 
*(/) = (/)> 

where (/) represents a net (u e ) such that u t = f for all e E (0, 1]. Then the map (13.121) 
is a linear embedding satisfying 

Ac^(Q) = cr- (3.14) 

Proof. This is a special case of the embedding constructed in [6] in the event that is 
bounded. □ 
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We will require Theorem 13.71 when we discuss how to solve (11.11) in Q(Cl). The ap- 
proach will be to reformulate the problem ( ll.lt as a differential equation in Q(Q) by 
utilizing a slight variation of the above embedding. But before we discuss how to do 
this, we must first define what we mean by a differential equation in Q(£i). 

3.4. Nets of Semilinear Differential Operators. We begin by defining a semilinear 
differential operator on £7(0). Our construction strongly resembles the construction by 
Mitrovic and Pilipovic in lfT2ll . For e < 1, if (a l e j ), (b\) E £ M (Q), we obtain a net of 
operators by defining A e to be 

K K 

A t u e = -DiiagDju) + J2 b > ni = —ai j DiDjU e - (D ia i j )(D jUe ) +J2 b >e) ni , 

i i=l 

where rii E Z. Under certain conditions, we can view a net of operators of the above 
form as an operator on Q(Cl). Here we determine these conditions, which will guarantee 
that this net of operators is a well-defined operator on 

Given an element u in Q{Vt), we first need to ensure that (A e u e ) E £m{®)- Based on 
how derivations and multiplication are defined in Q{Vt), the only serious obstacle to this 
is if rii < for some i < K. Therefore, we must guarantee that the element {{u e ) TH ) is a 
well-defined representative in Q{Vt) if rii < 0. It suffices to ensure that u = [(u e )] has an 
inverse in Q{Q). This is true if for each representative (u e ) of u, there exists e E (0, 1] 
and m E N such that for all e E (0, e ), hii^ |« e (^)| > Ce m . See (6l for more details. 
So u E G(Vt) must possess this property in order for the above operator to have any 
chance of being well-defined. For the rest of this section we assume that u satisfies this 
condition. 

Now suppose (a^J), (b e ) in £ M {Vt), and let 

N K K 

A e u = -Y, Dift'Dju) + Y^tu ni = K j DiD jUe - {D^i){DjU e ) + ^V e {u e ) n \ 

i,j=l i i=l 

We say that (A e ) ~ (A e ) if (a*-» - a^ j ), (6* - V e ) E A/" S (H). Then (A e ) ~ (A e ) if and only 
if (A e u e — A e u e ) E N(Vt) for all (u e ) E £m{&) due to the fact that the above operators 
are linear in (a l e j ) and (b\). 

Let A be the family of nets of differential operators of the above form and define 
Aq = A/ ~. Then for A E Aq and u E £m{&), define 

A : g(H) g(U) by Au = [A e u e ], 

where 

K 

[A e u e ] = [-0%] [DiDjUe] + [-D ia i j ][D jUe ] + Y} h % u Tl ( 3 - 15 ) 

i=i 

Using this definition, A E Aq is a well-defined operator on Q(Vl). We summarize this 
statement in the following proposition. 

Proposition 3.8. Aq is a well-defined class of differential operators from Q{VL) to Q{Q). 

Proof. Based on the construction of Aq, it is clear that for a given representative (u e ) of 
u E G(£l), (A e u e ) and (A e u e ) represent the same element in Q(Q). Furthermore, given a 
representative (A e ) of A , we also have that [A e u e ] = [A e u € ] for any two representatives 
of u E G(Vt). To see this, we first observe that for each e, every term in A e u e is linear 
except for the (u e ) Ui terms. So to verify the previous statement it suffices to show that for 
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each £ Z, ((w e ) ni ) = ((w e ) ni ) + {%), where (rj e )_e Af(tt). Given that [(u e )} = [(u e )} 
in G(Q), we have (u e ) = (u e ) + (r] e ) for (r] e ) G J\f(Q). For fixed e, n, G Z + , 



j=0 

where rj e consists of the summands that each contain some nonzero power of i] e . Clearly 
the net (rj e ) G jV(O). If rij G Z~, then for a fixed e, 

1 1 1 



for 



K + 77 e ) |nil x)^Q ( l 7 l )K) i (^) |nil " j K) |ni1 

By looking at the difference 

(u ) n > - = ^ = fi 

we see that the net ((u € ) ni ) = ((u € ) ni ) + (rj e ), where (f) € ) G Af(£l). Therefore for any 
u G G(Vl) possessing an inverse, and any A G Aq, the expression Au = [A e u € ] G Q{Vl) 
is well-defined. □ 

3.5. The Dirichlet Problem in Using the above definition of A, we can now 

define our semilinear Dirichlet problem on Q{Vt). Let u, p G Q{Vt) where O C M" is 
open, bounded and of C°° -class. Then let E be a total extension operator of Vt such that 
for / G C°°(fi), £/ G C°°(M n ) and = / See ffl for details. Using £ we may 
may define u\qq = p\gn for elements u, p G </(fi) if there are representatives (u e ) and 
(p e ) such that 

Ueldn = Pe\an + ^elan, 
where n e is a net of C°° functions defined in a neighborhood of d£l such that 

sup \n e (x)\ = o{e a ) Va G E. (3.16) 

x£dCl 

This will ensure that u\gn = p\an does not depend on representatives |[T2|. From now 
on we will abbreviate u\dn = p\gn by = P- With this definition of boundary 
equivalence, for a given operator A e Ao, the Dirichlet problem 

Au = inQ, (3.17) 

«|an = P, 

is well-defined in Q{Q). Now we state the conditions under which the above problem 
can be solved in Q{Vt). 

4. Overview of the Main Results 

We begin this section by stating the main existence result for the Dirichlet problem 
(13.171) . Let A G Aq be an operator on Q(Vt) defined by (13.151) . Also assume that the coef- 
ficients of A have representatives (a l J ) , (b\) G Em (ty that satisfy the following properties 
fore G (0,1): 

o« = af , al'Uj > A e |e| 2 > C ie a |e| 2 , (4-1) 
for each k eN, |a* J Ifc+i^fi, |&e|fc,a : n < A M < CV 6 , 
6* < -C 3 e c , {nj : n» < 0} 0, m = min{ni : n, < 0} 
&f > C*4e d , {rij : n« > 0} 7^ n^; = max{nj : 7ij > 0}, 
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where Ci,C 2 ,C 3 and C 4 are positive constants independent of e. Then the following 
Dirichlet problem has a solution in Q(Q): 

Au = [A £ u e ] = in ft, (4.2) 

u\an = P- 

We summarize this result in the following theorem, which will be the focus of the re- 
mainder of the paper: 

Theorem 4.1. Suppose that A : Q{Q) — > Q{Vt) is in Aq and that the conditions of (14.11) 
hold. Assume that for each e £ (0,1] that for some 1 < % < K, b\ is non-constant. 
Furthermore, assume that p £ £?(ft) has a representative (p e ) such that for e < 1, 
Pe > Ce a for some C > and a £ R. Then there exists a solution to the Dirichlet 
problem (|4T2l in £(ft). 

Proof. The proof will be given in Section |6j □ 

Remark 4.2. We can actually weaken the assumptions in (14.11) so that the conditions on 
the representatives (a l J), (&f ), (p e ) only have to hold for all e £ (0, eo) for some 
eo £ (0, 1). Suppose that this is the case, and that using these conditions we are able to 
show that for all e £ (0, eo), there exists u e that solves 

A e u e = in ft, (4.3) 

Ue\dU = Pe- 
lf U e satisfies the additional property that for all k £ N, there exists some e £ (0, eo), 
C > 0, and a £ R such that for all e £ (0, e ), |w e |fc, a < Ce a , then we can form a solution 
(v t ) £ £j\,f(ft) to (14.21) by defining v e = u t for e £ (0, eo) and v e = u eo for e £ [eo, 1]. The 
solution theory that we develop to prove Theorem \4. l\ with the stronger conditions (14.11 ) 
will also imply the existence of the partial net (u e ) of solutions to (14.31 ) in the event that 
the constraints outlined in (14.11 ) only hold for e £ (0, eo) C (0, 1). We will require this 
fact when we consider how to embed and solve (ll.lt in Q{Vt) later on in Section \43\ 

We begin assembling the tools we will need to prove Theorem 14. 1 L The first tool we 
need is a method capable of solving a large class of semilinear problems. The method of 
sub- and super-solutions meets this need, and we discuss this process of solving elliptic, 
semilinear problems in the following section. 

4. 1 . The Method of Sub- and Super-Solutions. In Theorem |431 below, we state a fixed 
point result that will be essential in proving Theorem |4.1[ This fixed point result is known 
as the method of sub- and super-solutions due to the fact that for a given operator A, the 
method relies on finding a sub-solution u_ and super- solution u + such that u_ < u + . 
A large part of this paper is devoted to finding a net of positive sub- and super-solutions 
for (14.21) and establishing growth conditions for them. In the proof below, let 

Lu = -Di(a ij Dju) + cu, (4.4) 

be an elliptic operator where 

a« = a 3i , a 'U, > A|£| 2 and a l \ceC°° (ft) . 

We now state and prove the sub- and super-solution fixed point result for these assump- 
tions. 

Theorem 4.3. Suppose ft C R n is a C°° domain and assume f : ft x R + — > R is in 
C°°(ft x R+) and p £ C°°(ft). Let L be of the form (l44b . Suppose that there exist 
functions U- : ft — > R and u + : ft — > R such that the following hold: 
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(1) u„,u + eC°°(TL), 

(2) < u-(x) < u+(x) Wx e n, 

(3) Lu- < f(x,u-), 

(4) Lu + > f (x,u+), 

(5) U- < p on dVL, 

(6) u + > p on dVt. 

Then there exists a solution u to 

Lu — f(x,u) onfl, (4.5) 
u\dn = P, 

such that 

(i) u e C°°(U), 

(ii) u-(x) < u(x) < u+(x). 

Proof. The general approach of the proof will be to construct a monotone sequence {u n } 
that is point-wise bounded above and below by our super- and sub-solutions, u + and 
U-. We will then apply elliptic regularity estimates and the Arzela-Ascoli Theorem to 
conclude that the sequence {u n } has a C°°(fi) limit u that is a solution to 

Lu = f(x,u) on f2, (4.6) 

u\on = P- 

Given that u_(x), u + (x) e C°°(Q), the interval [min U-{x), max + u + (x)\ C R + is 
well-defined. We then restrict the domain of the function / to the compact set K = 
n x [mmu-(x), max + u + (x)\. Given that / e C°°(£l x R+), it is clearly in C°°(Vl x 
[mmu-(x), max + u+(x)]) and so the function | d ^'^ \ is continuous and attains a max- 
imum on K. Denoting this maximum value by m, let M = max{m, — inf x6 ^ c(x)}. 
Then consider the operator 

Au = Lu + Mu, 

and the function 

F(x,t) = Mu + f(x,t). 

Note that this choice of M ensures that F(x, t) is an increasing function in t on K and 
that A is an invertible operator. Also, we clearly have the following: 

A(u) = F(x, u) <^ Lu = f(x, u), (4.7) 
A{u_) < F(x,uJ) L(u_) < f(x,u_), (4.8) 
A(u + ) > F(x, u + ) ^ L(u + ) > f(x, u + ). (4.9) 

The first step in the proof is to construct the sequence {u n } iteratively. Let u\ satisfy the 
equation 

A(ui) = F{x,u-) onQ, (4.10) 
ui\an = P- 

We observe that for u, v G Hq(Q), the operator A satisfies 

c i\\ u \\m(n) < ( Au , u ), and (Au,v) < ||w|||ri(n)lklllri(n)> 

where 

(u,v) = / uvdx, and (Lu,v) = / {a 13 DjuD(v + cuv)dx. 
Jn Jn 
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Therefore the Lax-Milgram theorem implies that there exists a weak solution U\ G H 1 (ft) 
satisfying U\ — p G Hq(Q). Given our assumptions on F(x, t) and p, F(x, u + ) G H m (il) 
and p G H m (VL) for all m G N. Therefore, by standard elliptic regularity arguments, 
u\ G H m (Q) for all m G N. This, the assumption that ft is of C°°-class and the assump- 
tion that a lj , c, p G C°°(ft) imply that «i G C°°(ft) and «i = p on <9ft. Therefore, we 
may iteratively define the sequence {«,} C C°°(ft) where 

A{uj) = F(x,uj-i) onft, (4.11) 
Uj\dn = P- 

The next step is to verify that the sequence {uj} is a monotonic increasing sequence 
satisfying u- < U\ < • • ■ < %_i < Uj < • - • < u+. We prove this by induction. First 
we observe that 

A(u^ -ui) < F(x,u-) -F(x,u-) = on ft, (4.12) 

(u_ - Mi) |an < 0. 

Therefore, by the weak maximum principle, U- < U\ on Vt. Now suppose that Uj-\ < Uj. 
Then 

A(uj — Uj+i) = F(x, Uj-i) — F(x, Uj) < on ft, (4.13) 

(wj - u j+1 )\asi = 0. 

given that F(x, t) is an increasing function in the variable t and Uj-i < Uj. The weak 
maximum principle again implies that Uj < Uj + i, so by induction we have that {uj} 
is monotonic increasing sequence that is point- wise bounded below by w_(x). Now we 
show that our increasing sequence is point-wise bounded above by u + (x) by proceeding 
in a similar manner. Given that w_ < u + and u + is a super-solution, we have that 

A(ux -u+) < F(x,u_) -F{x,u + ) < on ft, (4.14) 

(wi - u + )\ d n < 0. 

The weak maximum principle implies that u x < u + . Now assume that Uj < u + . Then 

A(u j+1 - u+) < F(x, uj) - F(x, u + ) < on ft, (4.15) 

(u j+ i - u + )\ an < 0, 

given that F(x, t) is an increasing function and Uj < u + . So by induction the sequence 
{uj} is a monotonic increasing sequence that is point- wise bounded above by u+(x) and 
point- wise bounded below by u-(x). 

Up to this point, we have constructed a monotonic increasing sequence {uj} C C°° (ft) 
such that for each j, Uj satisfies the Dirichlet problem (14.1 II) and is point- wise bounded 
below by u_ and above by u + . The next step will be to apply the Arzela-Ascoli theorem 
and a bootstrapping argument to conclude that this sequence converges to u G C°°(ft). 
We first show that it converges to u G C (ft) by an application of the Arzela-Ascoli The- 
orem. Clearly the family of functions {«, } is point- wise bounded, so it is only necessary 
to establish the equicontinuity of the sequence. Given that each function Uj solves the 
problem (14.1 II) . by standard L p elliptic regularity estimates (cf. [|5]|) we have that 

||«il|w a .f < C(\\Uj\\ L v + \\F(x, %_i)||lp)- 

The regularity of F(x,t) and the sequence {uj} along with the above estimate and 
the compactness of ft x [inf u_, sup w + ] imply that there exists a constant N such that 



18 



M. HOLST AND C. MEIER 



\\F(x,Uj_i)\\ LP < N for all j. Therefore, if p > 3, the above bound and the fact that 
w_ < Uj < u + imply that for each j G N, 

|wj|i ja; n < C||w||iy2, P < oo, 

where a = 1 — -. This implies that the sequence {uj} is equicontinuous. The Arzela- 

Ascoli Theorem then implies that there exists a u G C(fi) and a subsequence {uj k } such 
that M Jfe — > u uniformly. Furthermore, due to the fact that the sequence {uj} is monotonic 
increasing, we actually have that Uj — > u uniformly on Q. Once we have that Uj — > it in 
C(fi), we apply L p regularity theory again to conclude that 

\uj - u k \i, a -n <C||uj - u k \\ W 2, P (4.16) 
<C'(\\uj - u k \\ LP + ||F(a?,«j-_i) - 

Note that the above estimate follows from the fact that Uj k+ \ — Uj l+ \ satisfies 

A(uj — life) = F(x, Uj-i) — F(x, Uk-i) on fi, (4.17) 
(uj - u k )\dn = 0. 

Given that Uj — > u in C(Q), (14.161) implies that the sequence {uj} is a Cauchy sequence 
in C 1 (f2). The completeness of C 1 (fi) then implies that this subsequence has a limit 
v G C 1 (fi), and given that Uj — > u in C(f2), it follows that u = t>. Similarly, by 
repeating the above argument and using higher order L p estimates we have that 

\uj - u k \ 2t a-n <C{\\uj - u k \\ W 3, P ) (4.18) 
<C'(\\uj - u k \\ w i, P + \\F(x,Uj-i) ~ F(a?,« fc _i)|| W ri,p), 

where — >■ u in C 1 (r2) as fc — >■ oo. Again, (14.181) . the regularity of F and the fact that 
Uj — y u in C 1 (fi) imply that the sequence {uj} is Cauchy in C 2 (Q). A simple induction 
argument then shows that u G C°°(fi). 

The final step of the proof is to show that u is an actual solution to the problem (I4.5I ). 
It suffices to show that u is a weak solution to the above problem. It is clear that u = 
p on dVt, so we only need to show that u satisfies (14.51) on Vt. Fix v G Then 
based on the definition of the sequence {uj}, we have 



(a lJ DjUjDiV + Mujv)dx = / (/(x, m 3 _x) + Muj_i)vdx. 
i Jn 

As Wj — )■ w uniformly in C(Q), we let let j — > oo to conclude that 

(a lj DjuDiV + Muv)dx = I (f(x, u) + Mu)vdx. 



Upon canceling the term involving M from both sides, we find that u is a weak solution. 

□ 



4.2. Outline of the Proof of Theorem 14.11 Now that the sub- and super- solution fixed 
point theorem is in place, we give an outline for how to prove Theorem 14. II 

Step 1 : Formulation of the problem. We phrase (|4.2I) in a way that allows us to solve 
a net of semilinear elliptic problems. We assume that the coefficients of A and 
boundary data p have representatives (a l J), (b\), and (p € ) in £m(Q) satisfying the 
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assumptions (I4.ll ). Then for this particular choice of representatives, we solve 
the family of problems: 

N N 

A e u e = -J2 DiK j D jUe ) + tie u 7 = in fi > ( 4 - 19 ) 

i,j=l i 

u € \an = Pe- 

Then we must ensure that the net of solutions (u e ) E £m(Q) and ensure that (14.191 ) 
is satisfied for other representatives of A, p, u. 
Step 2: Determine L°° -estimates and a net of generalized constant sub-solutions and 
super- solutions. We determine constant, a priori L°° bounds such that for a pos- 
itive net of solutions (u e ) of the semilinear problem (14.191 ), there exist constants 
oi, 02 G R, C\, C2 > independent of e G (0, 1) such that 

Cie ai < a £ < u e < p £ < C 2 e a2 . 

These estimates are constructed in such a way that for each e, the pair a € , [3 e are 
sub- and super-solutions for (14.191) . 
Step 3: Apply fixed point theorem to solve each semilinear problem in (14. 191 ). Using the 
sub- and super-solutions a e , f3 e , we apply Theorem 14. 31 to obtain a net of solutions 

(u e ) G C°°(Jl). 

Step 4: Verify that the net of solutions (u e ) G £m(Q). Here we show that the net of so- 
lutions satisfies the necessary growth conditions in e using the growth conditions 
on the sub- and super- solutions and Theorem l3.ll 

Step 5: Verify that the solution is well-defined. Once we've determined that the net of 
solutions (u e ) G £m{&), we conclude that [(u e )] G G(Cl) is a solution to the 
Dirichlet problem (14.21) by showing that the solution is independent of the repre- 
sentatives chosen. Note that most of the work for this step was done in Proposi- 
tion m 

We shall carry out the above steps in our proof of Theorem 14. II in Section |6] We still 
need to determine a net of sub- and super- solutions for (|4.1I) . which we do in Section |5l 
But before we move on to this and the other steps in the above outline, we briefly return to 
the motivating problem (ll.lt by discussing how to embed a problem with distributional 
data into Q{Vt). 

4.3. Embedding a Semilinear Elliptic PDE with Distributional Data into Q (Q). Now 
that we have defined what it means to solve a differential equation in Q(Q), we are ready 
to return to the problem discussed at the beginning of the paper. We are interested in 
solving an elliptic, semilinear Dirichlet problem of the form 

N K 

- Di{a ij D jU ) + = infl, (4.20) 

i,j=X i=l 

u\an = p, 

where a 1 - 7 , b l and p are potentially distributional and n, G Z for each i. If we can formu- 
late this problem as a family of equations similar to (14-. 1 9b . then it can readily be solved 
in Q(tt) by Theorem 14.11 What we require is a method of embedding our problem into 
Q(Cl). The key to constructing such an embedding will be Theorem 13 .71 

For the following discussion, let CC £1' C M. n be open and bounded. Additionally 
assume that is of C°°-class. Given a distribution u E V'(Cl'), define the restriction u\a 
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to be u restricted to test functions p G T>(Q). Let V'(Q')\q C V(Q) be the space of 
distributions obtained in this way. Additionally, given an element v G Q s (f2'), define 

v\n = (ve\n)+^m- 
Then we have the following proposition: 

Proposition 4.4. Let i denote the embedding defined in Theorem \3. 71 Define the follow- 
ing map 

a :C°°(Q) -^Q(Q), (4.21) 
*(/) = (/), 

where (/) w a ne? (u e ) smc/? that u e = f for all e G (0, 1]. TTzen 

i : X>'(fi')|n -> (/(H), (4.22) 
i(u\ u ) = (i(ti)lTy), 

w a linear embedding such thati\ Coo ^ = a. 

Proof. That « is linear and well-defined follows from the properties of restriction maps 
and the fact that i is linear and injective. To verify that i\ C oo/m = o, we use that fact that 
Q is of C°°-class, which implies that there exists a total extension operator E such that 

E : C°°{Q) -> C°°(M n ), (4.23) 
(Ef)\ u = f foreach/GC°°(fi). (4.24) 

See [1] for details. Then for each / G C°°^f), f = (Ef)\ n . Clearly Ef G V'{Q!), and 
therefore, C°°(VL) C The fact that «| COO (q) = o then follows from Theorem l3.7l 

and the definition of i. Finally, we need to show that i is injective. Suppose that i(u) = 
for some u G P'(O') |n. Then if (j> e and ^ are the same as in Theorem 13 .71 we have that 

(<(«)i7y) = (((e«)*0 e )yG^(s)- 

This implies that 

((^w) * e ) — >■ uniformly on VI. 
Therefore, for any ip G T>(Cl), we have that 

(u, if) = (u, £<p) = (£u, if) = lim ((£u) * <f> e , if) = 0. (4.25) 

£-)>0 

So i is injective. □ 

Now suppose that we are given a problem of the form (|4.20l) . where the terms a % \ b % 
and p are in ViVt'). Then we may use Proposition 14.41 to embed the coefficients a % \ b 1 
and p into Q . We will denote a representative of the image of each these terms in 
Q(Q) by (a l J), (b l e ) and (p e ). Then for a choice of representatives, we obtain a net of 
problems of the form (14. 191 ). 

In order to solve this net of problems by using Theorem 14. 1[ we need there to exist 
a choice of representatives (a*/), (&*) and (p e ) that satisfy the conditions specified in 
(I4.ll ). While these conditions might seem exacting, this solution framework still admits 
a wide range of interesting problems. This is evident when one considers the following 
proposition: 
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Proposition 4.5. Let rii E Z be a collection of integers for 1 < % < K. Assume that 
there exist 1 <i,j < K such that rii < and rij > and let 

n\ = min{nj : rii < 0}, and tlk = max{nj : rii > 0}. 

Suppose that a ij , b 1 , b K , p E C(tt') and b 2 ,--- , b K ~ 1 E V{9!), where Vt' C W 1 is an 
open and bounded set. Additionally assume that satisfies the symmetric, ellipticity 
condition and p > 0, b\ < and bx > in Q'. Then ifVL CC Q' is of C°° -class, the 
problem 

N K 

- Di{a ij Dju) + biuUi = in fi > ( 4 - 26 ) 

i,j=l i=l 

u\an = P, 

admits a solution in Q(Q). 

Proof. This follows from Proposition @31 Theorem 14. 11 Remark |4~21 and the fact that 
(a^*(p e ), (b 1 *^^, (b K *(f) e ) and (p*4> t ) converge uniformly to a JJ , b 1 , b K and p'mVl. For 
e sufficiently small, the corresponding problem (14.191) in Q (O) will satisfy the conditions 
specified in (14.11) . Therefore, Theorem 14. II and Remark l4~2l imply the result. □ 

With the issue of solving (14.201) at least partially resolved, we return to the task of 
proving Theorem 14. II We begin by establishing some a priori L°° -bounds for a solution 
to our semilinear problem (14.261) if the given data is smooth. 

5. Sub- and Super-Solution Construction and Estimates 

Given an operator A E A with coefficients satisfying (14.11) . our solution strategy for 
the Dirichlet problem (|4.2I) is to solve the family of problems (14.191) and then establish the 
necessary decay estimates. In order for this to be a viable strategy, we first need to show 
that (14.191 ) has a solution for each e G (0, 1). Given that rii < for some 1 < i < K, for 
each e, we must restrict the operator 

N K 
i,j=l i=l 

to a subset of functions in C°°(ri) to guarantee that A e is well-defined. In particular, for 
each e we consider functions u e G C°°(fi) such that < a t < u e < f3 t < oo for some 
choice of a e and (3 e . The first part of this section is dedicated to making judicious choices 
of a e and f3 e for each e such that a solution u t to (14.191 ) exists that satisfies a e < u t < f3 £ . 

Once a net of solutions (u e ) is determined, it is necessary to show that if (u e ) G £m(Q), 
then an operator A E A whose coefficients satisfy (14.11) is well-defined for (u e ). Recall 
that A is only a well defined operator for elements u E satisfying u t > Ce a for 

e G (0, e ) C (0, 1), a £ 1 and some constant C independent of e. This will require us 
to establish certain e-decay estimates on a e , which we do later in this section. 

5.1. L°° Bounds for the Semilinear Problem. We begin by determining the net of 
a priori bounds a e and (3 e described above. For now we disregard the e notation. In 
the following proposition we determine a priori estimates for a positive solution u to a 
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problem of the form 

N K 

- Di(a ij D jU ) + h%uni = in ft, (5.1) 

i,j i=l 

u\dn = p, 

where £1 C W 1 is connected, bounded, and of C°°-class, and a l \ b\ p G C(Q,) with p > 
in ft. 

Proposition 5.1. Suppose that the semilinear operator in (|5.1I) has the property that 
rii > for some 1 < % < K. Let nx be the largest positive exponent and suppose that 
b K (x) > in ft. Additionally, assume that one of the following two cases holds: 

(1) rii < for some l<i< K and if ^ = minjrij : rij < 0} , (5.2) 

then b± (x) < in fi. 

(2) uk is odd and < for all I < i < K. (5.3) 

Then if case (15.21 ) holds and u is a nonnegative solution to (15.11 ). there exist constants a 
and (3 such that < a < u < (3. If case (15.31 ) holds and u is a nonnegative solution to 
(15.11) . there exists a constant (3 such that < u < f3. 

Remark 5.2. Note that for the purposes of proving Theorem \4.1\ we are primarily con- 
cerned with case (15.21) . This is the case that we will focus on for the remainder of the 
paper. However, with a little extra work we could very easily generalize Theorem \4.1\ to 
allow for rii > Ofor all 1 < i < K and n K > odd. Then we could use case (15.31 ) to 
establish the necessary bounds. 

Proof. We first define a' and f3' in the following way. If case (15.21) holds, then let 

a' = sup { V sup bi(x)b n > < V6 G (0, c) 1 . (5.4) 

If case (15731) holds, let a' = 0. In either case ((5721 or ((573l . define 

(3' = inf | ^ inf bi(x)b ni > Wb G (c, oo) j . (5.5) 

Note that in all cases a' and (3' are always well-defined given the conditions on b 1 (x) and 
b K (x). Then define 

a = min{a', inf p(x)}, (5.6) 

(3 = m&x{(3', sup p(x)}. (5.7) 

xedn 

Based on these definitions of a and (3, if u is a nonnegative solution to (15.11 ) then it is 
easy to verify that the functions <\> = [u — (3) + and <p = [u — a) ~ are in if either 

(15721 or d57J) holds. Define the set 

y= {xeU\u>f3} 
if case (15721 or (15731 holds. If case (15721 holds, let 

y = {x e n I < u < a}, 

if case (15.31 ) holds, let 

y = {x G H I it < a}. 



GENERALIZED SOLUTIONS TO SEMILINEAR ELLIPTIC PDE 23 

Then if u G is a weak solution to (15.11 ), supp((/>) = y and supp(0) = y. We have 

the following string of inequalities for <p: 

C2\W HH n)<Ci\\V{{u-ar)\\l HQ) (5.8) 

< / a ij D j {{u-a)~)D j ({u-aY)dx 
Jn 

= / a l i Dj(u — a)Dj((u — a)~) dx 
Jn 

K 

= / (- y2bi(x)u n *)(u-a) dx < 0. 

J y i= i 

Similarly, we have the following string of inequalities for <p: 

C 2 ||^|||, 1{0) < C 1 ||V((u-/3) + )||| 2(0) (5.9) 

< / a^((M-/3) + )A((«-/3) + ) dx 



a 

a ij Dj(u- P)Di((u- dx 

n 

A" 

= /_(- y"h{x)u n *){u - ft) dx < 0. 

The above inequalities imply that if u is a positive, weak solution to the semilinear (|5.1I) . 
then w G [a, /3] where in case (15.21 ), cc > 0. □ 

Now that we've established L°°-bounds for solutions to (15.11) , we can apply these 
bounds for each fixed e to determine a net of bounds for the following net of problems: 

N K 

- D l al j D j u e + K U T = ° in n ( 5 - 10 ) 

i,j i=l 

Ue\dU = Pe, 

where (a l J), (b\), (p e ) G £m(Q) satisfy the following for all e < 1: 

a* = a*, agtti > Xe\e > C^ltf (5.11) 

\a*\ k , a ;n, \K\ k , a ;n<Ak,e<C 2 e a2 , Vk G N 

b\ < -C 3 e a; \ {rii : m < 0} ^ 0, m = min{nj : m < 0} 

bf > C 4 e a4 , {rij : m > 0} 7^ 0, n K = max{n^ : > 0} 

p e >C 5 e a5 , 

and Cx, • • • , C 5 are independent of e and a 1; • • • , a 5 G M. 

Proposition 5.3. Suppose that for each fixed e G (0, 1], w e z's a positive solution to (15.101) 
with coefficients satisfying (15.1 II) . 77zen there exist L°° -bounds a e and f3 e such that for 
each e, < a e < u e < (3 e . 

Proof. For each fixed e, if the assumptions in (15.1 II) hold, then case (|5.2I) of Propo- 
sition \5A\ is satisfied. Therefore, for each e G (0, 1], there exists a e and j3 e such that 

< a e < u e < (3 e . □ 
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5.2. Sub- and Super- Solutions. In the previous section we showed that if the data of 
(15.101) satisfies (15.111) and if u e E C°°(U) solves (15.101) for each e, then < a t < 
Ue < e - Now, for each e E (0, 1], we want to show that there actually exists a solution 
u e E C°°(f2) satisfying < a € < u e < (3 e . The key to proving this result lies in the fact 
that a t and f3 e are sub- and super- solutions to (15.101) for each e. 

Proposition 5.4. Let a t and (3 e be the bounds established in Proposition 15.31 Suppose 
that the coefficients in the net of problems (15.101) satisfy (|5.1 II) and that for each e some 
b\(x) is non-constant. Then there exists a net (u € ) E (C(f2)) 7 such that for each e, u e 
solves (15.101 ) and < a t < u e < (3 e . 

Proof To solve the above family of problems in (15.101) . we show that the net of L°°- 
bounds (a e ) and found in Proposition l5.3l is a net of sub and super-solutions to (15.101) . 
Then we verify that the interval [a<=,/3 e ] is nonempty and is a subset of R + . We will 
then be able to apply Theorem 14.31 to conclude that for each e, there exists a solution 

u e E C°°(Q). 

Fix e and let a' e and f3' e be defined by (|5.4I) and (|5.5I) respectively, and let 

a e = min{c^,inf p t (x)}, 

(3 t = max{/^, sup p t (x)}. 

By (1541) and the fact that p t > 0, we clearly have that a e > 0. Then (l5~4l) and the 
definition of a e imply that 

K K 

A e a e = Y,K{aeT l < J2 sn l b H a ^ ni ^ °> ( 5 - 12 ) 

i=l i=l 
&e < mf p e (x) < p e , 

which shows that a t is sub-solution for each e. Similarly, (15.51) and the definition of 
imply that 

K K 

(3 e > sup p € > p e , 
xedn 

which shows that (3 e is a super-solution for each e. 

Now that we've determined that the pair a e and /3 e are sub- and super-solutions of (15.101 ), 
we show that the interval [a e , f3 e ] is nonempty. Given the definition of a e and (3 e , it suf- 
fices to show that a' e < (3' € . We claim that a' € < (3' € provided that at least one b\(x) is 
nonzero and non-constant. This follows from the fact that if we define 

K 

7 £ = inf { V sup b i (x)d ni >0 We (c, oo)}, 

i=l xGil 

then we have that a' e < 7 e by (15.41 ) and (15.1 II ). Furthermore, for a fixed e, given the 
assumptions on 

K K 

i=i x£U i=i zefi 

But (15.51) and the above inequality clearly imply that 7 e < f3' e . Therefore a' e < f3' e and 
the interval [a e , f3 £ ] is a nonempty subset of M + . For each e E (0, 1], the hypotheses of 
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Theorem 14.31 are satisfied for the elliptic problem (15.1 It , so we may conclude that there 
exists a net of solutions (u e ) E (C 00 ^)) 1 that satisfy < a e < u € < f3 e for each fixed 
e. " □ 

The final task in this section is to show that an operator A E Aq, with coefficients 
satisfying (15.1 It , is a well-defined operator on any element u E £m(Q) satisfying 

ae<u t <(3 e Vee(0,l]. 

Recall that in Section 13.41 we determined that A is only well-defined for invertible u E 
Q(Cl). Therefore, it suffices to show that (a t ), (/3 e ) and (— ), (-3-) are generalized con- 
stants (13.31) . which we verify in the following lemma. 

Lemma 5.5. Let (a € ) and (/3 e ) be the net of sub- and super- solutions to (|5.10h determined 
in Section \5.1\ Suppose that the coefficients of (15. 101 ) satisfy (|5. 1 It . Then (a t ), (/3 e ),(^-), 

and are in C, the ring of generalized constants. 

Remark 5.6. Note that if (— ) E C, then this implies that there exists an eo E (0, 1), 
some constant C independent of e and a E R such that a e > Ce a foralle E (0, eo). Then 
if (u t ) E £m{Q) satisfies a t < u e < {3 e for each e, (^-) E C implies that u = [(u e )] is 
invertible in Q{VL). See Section U~4\ and ^ for more details. 

Proof. We need to show that there exists constants D\, D 2 independent of e and e E 
(0, 1) such that for all e E (0, e ), 

« £ > Dxe hl for some b x E R, 

fie < D 2 e b2 for some b 2 E R. 

So it is necessary to verify that there exists constants D\ and D 2 so that for e sufficiently 
small 

a' > LV\ and inf p e > D^ 1 , 

Pi < D 2 e h \ and sup p e < D 2 e b2 . 
xedn 

Given that (p e ) E £ M (P), 

sup p e < sup p e = 0{e b ) 1 

for some b E R. This and the assumption on (p e ) in (|5.1 It imply that we only need to 
obtain the necessary e-bounds on a' e and f3' e . 

For now, drop the e notation and consider a' defined in (15 .4t . For a given function /, 
define 

7 = sup{/(6)<0 V6G(0,c)}. 

1 cGM+ 

Given that 

a' = sup \ V sup bi(x)b nt < Mb E (0, c) I , 
it is clear that for another function f(y) such that 

K 

f(y) >J2 su ^ bi ^y ni on (°' c )' 

1=1 

if 7^ is defined and a' E (0, c), it must hold that 7 < a'. Let Ci = \{ni : > 0}| 
and C2 = : rii < 0}| and if C2 > 1, let n i2 = min{nj : ri\ < rii < 0}. Note that 
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Ci,C 2 > 1 based on the assumptions in (15. 1 II ). Then recalling that bi(x) < 0,bx(x) > 
correspond to the coefficients of the terms with the smallest negative and largest positive 
exponent of £\ bi(x)u n \ if sup^Q \bi(x)\ < A for each i, the following must hold for 

ye (o,i): 

K 

J2 su PH%)y ni < supfe^x)?/™ 1 + dA + (C 2 - l)Ay n *2. (5.14) 

Define 

V 2(C 2 -1)A 

if Ci > 1 and let d = 1 if C 2 = 1. Then let c = min{l, c?}. The definition of c implies 
that 

(C 2 - 1)A^ < - SUp ^ 6l(x) y"% 
for all 1/ G (0, c). So for y G (0, c), 

fVpfc^)^ < SUP ' e f (g) ^ + dA = /(„). 

Then if a' G (0, c), a' > 7 . Given that f(y) is a monotone increasing function on R + , 
7^ is the lone positive root of f(y). Thus, 



7, 



1 

-"•1 



/ V 2dA 
which implies that if a' G (0, c), 

2d A 

Similarly, for a fixed e G (0, 1), define 



a' > 



(L 



swp x& u,{b\(x))\^2 



H 

J 



2(C 2 - l)A e 

if C 2 > 1 and let d e = 1 if C 2 = 1. Let c e = min{l, d e }. Then for y G (0, c e ), we have 
that 

(C 2 - l)A e y^ < _!^§|5M y m 

So the above arguments imply that if a' e G (0, c e ), then > 7 and 

1 



Given the assumptions on b\(x) and A e in (15.1 II) . in this case we have that a' e > Ce a for 
some constant C > 0, a G R and e sufficiently small. Now we must show that c e > Ce a 
for some constant C > 0, a G R and e sufficiently small in the event that a£ ^ (0, c e ). 
It suffices to show that d e > Ce a in the event that C 2 > 1. But clearly, for e sufficiently 
small 



e ~ ' "2(d-l)A e 1 " ( • 
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given the assumptions on b\ and A e in (15.1 II ). Therefore a' e > D\e a for some constant 
Di > 0, a G R and e sufficiently small. 

Now we determine bounds on the net (j3' e ). Again, we temporarily drop the e and only 
consider f3' . Recall that 

13' = inf inf bi{x)b ni > Mb G (c, oo) 

For a given function /(y), define 

j f = inf {/(6) > V6 G (c,oo)}. 

Then if /(y) < J2f=i su Pzef2 bi(x)y ni on some interval (c, oo) and j3' G (c, oo), it must 
hold that 7 j > if 7^ is defined. Let C\, C 2 be as before and let = max{nj : < 
7ii < n K } if Cx > 1. If y > 1, then 

if 

V inf 6i(x)?/ ni > inf (^(a^y™* - (Ci - l)Ay ni i - C 2 A. 

Now define 

2(d - 1)A 



d 



if Ci > 1 and let c? = 1 if C± — 1. Let c = max{l, c?}. Then our choice of d ensures that 
if Ci > 1, then 

-( Cl - l)Ay-n > _^(M*))y~ 
and that for y G (c, 00), 



K 



J] supfc^y- > ^faW yn* - C7 2 A = /(y). 



i=l 



2 



So if 0' G (c, 00), (3' < 7j, where 7^ is the lone positive root of / on R + given that / is 
monotone increasing on this interval. So if /3' G (c, 00), 



0'<l f 

By defining 



2C 2 A 



i^enCM^)) 



d e = [ ~~^~n~F7~\~ I 1 anc ^ c e = max{l, e? e } (5.15) 



inf bf-(x) 

and applying the above argument for (3' to the net ((3' e ) for each fixed e, it is clear that if 

f3' e G (c e , 00), then 

/ 2C7 2 A £ \ 4 

given the assumptions on 6^ and A e in (|5.1 II) . 

Now assume that f3' t (c e , 00). Then it suffices to show that if C\ > 1, then d e < Ce a 
for e sufficiently small, some positive constant C and a G R. But again, this is clearly 
true given the assumptions (|5.1 II) and the fact that 



2(Ci - l)A t 
inf bf (x) 
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□ 

6. Proof of the Main Results 

We now prove Theorem 14 . 1 1 using the results from Section [5l For clarity, we break the 
proof up into the steps outlined in Section l4~2l 

6.1. Proof of Theorem HI 

Proof. Step 1: Formulation of the problem. For convenience, we restate the problem and 
the formulation that we will use to find a solution. Given an operator A E Aq, 
defined by (13.151) , we want to solve the following Dirichlet problem in 

Au = infi, (6.1) 
u\an = p . (6.2) 

We phrase (16.11) in a way that allows us to solve a net of semilinear elliptic 
problems. We assume that the coefficients of A and boundary data p have repre- 
sentatives (a\ j ), (&*), and (p e ) in £m(&) satisfying the assumptions (14.11) . Then 
for this particular choice of representatives, our strategy for solving (16.11) is to 
solve the family of problems 

N N 

A e u e = -J2 D^DjUe) + K U T = in ft, (6.3) 

i,j=l i 

u e \dn = Pe, (6.4) 

and then show that the net of solutions [u t ) E £ M (fi). 

Step 2: Determine L°° -estimates and a net of sub-solutions and super-solutions. In Sec- 
tion |51 we concluded that for each e, the pair a e and f3 e determine sub- and super- 
solutions to (|6.3I) such that < a € < (3 e . Furthermore, in Lemma 1531 we con- 
cluded that there exist C\, C 2 > and a 1; a 2 E K. such that for e sufficiently 
small, the nets (a e ) and (f3 e ) satisfy C\e ai < a e < (3 e < C 2 e a2 , thereby verifying 
that (a e ), ((3 e ), (^-), (j-) E C, the ring of generalized constants. 

Step 3: Apply fixed point theorem to solve each semilinear problem in (14.191) . This fol- 
lows from Proposition 15 .41 We briefly reiterate the proof here. We simply verify 
the hypotheses of Theorem |4.31 For each fixed e we have sub- and super-solutions 
a € and f3 e satisfying < a € < f3 e and a 1 /, b\, p e E C 00 (f2) satisfying (15.111) . Fi- 
nally, f2 is of C°°-class and the function 

K 

f(x,y) = -J2K(x)v ni e C°°(QxR + ), 

8=1 

so we may apply Theorem 14 .3 1 to conclude that there exists a net of solutions (u e ) 
to (15.101 ) satisfying < a e < u e < f3 e . 
Step 4: Verify that the net of solutions (u e ) E £m(Q)- Now that it is clear that a solution 
exists for (15.101 ) for each e E (0, 1], it is necessary to establish estimates that 
show that the net of solutions (u e ) is in That is, we want to show that for 

each k EN and all multi-indices \(3\ < k, there exists a6l such that 

sup{\DP Ue (x)\} = 0{e a ). 
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By standard interpolation inequalities, it suffices to show that for 7 6 (0, 1) and 
each k G N, there exists anaGl such that 

kl*, 7; n = 0{e a ). 

By Theorem [3TT1 we have that if u € is a solution to (15.101) with coefficients satis- 
fying (15.111) . then 

\u t \2, r & 

- C {f) (l^loin + l^l^n + zZl 6 ^)"*! ^ )- (6 - 5) 

Observe that 

K'|o, 7; n < K ! |o ; f7 +^[«e]o,7;n|We|^n 1 (6.6) 



if rii > and 

|< l |o, 7; n < |< s |o;Q + | z^Ta (-^)Ho, 7 ;nl^l^n" 1 » ( 6 -?) 



l Ue l0;O 



if rij < 0. The above inequality implies that 



|«e|2,7;fi 

A 

+ 22 \K{ x )\o, T ,n( c i( a ei Pei n i) + c 2(rii, a e , (3 e )\u e \ 0tr> n)), 

i=l 

where 

Ci(7ii,a e ,^ e ) = and C 2 (ni,a e ,/3 e ) = if n f > and 

Ci(7ii, a e , = a" 1 and C 2 {n h a e , e ) = if ti< < 0. 

Application of the interpolation inequality 

I W e 1 0,7 < C(8~ X \u e \ Q + 5 e \u e \ 2 ,j) } 

where 5 e is arbitrarily small and C is independent of 5 £ , implies that 

A^ 3 



\ue\2, r ,n < C ( — ) (|w e | 0; n + \Pe\2, r ,n (6.9) 



A 



+ l &l e( x )lo,7;n(Ci(7ii,a e ,/3 e ) 



i=i 



+ C 2 (n i ,a e ,/3 e )(C(5 e 1 |« e | 0; n + tfekla^n))))- 

Therefore, 

1 -<J e (^j ^\K( X )\o^ C 2( n i, a e,Pe)j Kh^O (6.10) 



i=l 

3 



< C ^ ) (|tt e | 0; n + |P e |2, 7 ;Q 
A' 



+ ^ |& l e (x)|o, 7 ;f7(C , i(n i ,a e ,/3 e ) + C 2 (nj, a e , /3 e )6 e 1 \u e \ 0; n))- 



8=1 
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But given the assumptions on A e , A e , the bounds previously established for the 
nets (a £ ) and (f3 e ) in Lemma 1531 and given that (b\(x)) G Em{&), there exists 
e G (0, 1), a G M and C > such that for all e G (0, e ), 



'A 1 '• 



Therefore, choosing 



Y,\ h \{ x )Wi C ^^^e)<C^. 



8=1 



6. 1 



2Ce a ' 

it is clear that for e G (0, e ), 

3 



Ue\2, r ,n < C (J^j (\u e \o;n + \Pe\2, r ,n (6.11) 

K 

+ £ \K( x )\o, r ,n(Ci(ni, a e , (3 e ) + C 2 (n il a £ , e , e a )\u e \o-,n)) ■ 



8=1 



Given that (a € ), (/3 e ) G C, a e < u e < /3 e and (p e ), (b\) G £^-($1), the above 
inequality implies that for some o£l, 

K\ 2 , r ,n = 0{t a ). 

Now we need to utilize the e-growth conditions on |u e |2, 7; c and induction to show 
that for any k > 2 that 

\ue\k, r ,n = 0{t a ) for some aeR. (6.12) 

Let (u e ) be a smooth net of solutions to (|6.3I) and additionally assume that (16.121) 
holds for all j < k. Let v be a multi-index of length k — 1 . Then by differentiating 
both sides of (16.31) . we see that for each e, w e satisfies the Dirichlet problem 

N K 

J2 D^-D^D^)) = -^D v {V € u?) infi (6.13) 

i,j=l i=l 

D u u e \ 9 Q = D u p e . 

Rearranging the above equation and applying the multi-index product rule we 
find that 

N N 

£ a^'AiP^e) = ~J2 D "(( D i a ?)( D 3 u t)) (6.14) 

i,j=l i,j=l 

N 



E E £sWW D 

r + ^i— 
(7 7^ 



— ' cnu 

i,j = l (y+ti=v 



K 



(TILL 

1=1 a+^i=u 
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Therefore, we may apply Theorem [3J] to (16.141) to conclude that for an arbitrary 
multi-index v such that \v\ = k — 1, 

\D u u e \ 2il , n <c(^ (\D v u e \ .,n + \D v Pe \ 2iT , Q (6.15) 

N 

+ \J2 DV (iD i a^)(D j u e ))\ 0jJ . Q 

N , 
+ EE ^|^X J |0, 7 ;d^A^|0, 7 ;n 



K 

V 



+ EE -y-Tl^lo,7;d^(K) n! )lo, 7; n) 



(Till 

1=1 (T+fJ,=V 



<C\f) {\D v u e \ Q , n + \D» Pe \ 2>7 . Q 



N i 



cnu 

l,J = l (T+(J,=U 

N ! 



i,j = l <T-\-fi = U 

K , 

+ E E -r3p^lo,7;np CT ((« e ) ni )lo, 7; n). 
* — ' * — ' o"!u! 

1=1 CT+/^=^ 

By our inductive hypothesis and the assumptions on the coefficients, it is imme- 
diate that every term in the above expression is 0(e a ) for some oGl except for 
the last term. So to show 

\D u u e \ 2 , r ,Q = 0(e a ) for some aGl, 

it suffices to show that 

E E ^D^KW r ,n\D^{u e )^)\^ = 0{e a ) for some a G R. 

i=l <r+fj,=i; 

Given that b\ e S M (Jl) for each l<i<K, 

|£W e |o, 7; n = C>(e a ) for some aGl 

Therefore, it is really only necessary to show that for any multi-index a, such that 
Wl = j < k — 1, that there exists ana G M such that 

\D°((u e D\ , r ,n = O(e a ). 

But observe that D a ({u e ) ni ) is a sum of terms of the form 

(u e ) ni - m D ai u e D a2 u e ■ ■ ■ D a ™u e , 
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where a 1 + a 2 + ■ ■ ■ o m = a and m < j < k — 1. This follows immediately from 
the chain rule. Therefore we have the following bound: 



iZT^nio,^ < (nO|(iO n< -\ 7 ;d£ a «e|o, 7 ;n (6.16) 

(T\+(T2=cr 

■ \D^u e \ , r M Da2u Xr,n + --- 

+ J2 — j — r — r(Ofai-i) 

o\\a-2}. • • • a A 

fl+(T2H 1-^=°" 

• • • (rii - j)\(u e ) n >- j \o,r,n\D ai u e \ ,r,n 

■ ■ ■ \D aj u e \ , r ,Q. 

Using (16.61) and (16.71) . for each m < j we may bound the terms of the form 
|(w e ) ni_m |o, 7; Q using |ii e |o, 7; n, cn' t and f3' e . Then our inductive hypothesis and the 
growth conditions on (a' e ) and (f3' e ) imply that 

\D a ((u € ) ni )\o,rV = °^ a ) for some a e K 
This implies that 

|-D"w e |2, 7; n = C(e a ) for some o6l. 

As was an arbitrary multi-index such that \u\ = k — 1, this implies there exists 
fl£M such that 

|«eU+i, 7; n = 0(e a ). 

Therefore, (u e ) E £ M (Q)- 
Step 5: Verify that the solution is well-defined. Proposition 13 .81 and the definition of the 
Dirichlet problem in Q(Vl) given in Section [331 imply that [(u e )} is indeed a solu- 
tion to the problem 

Au = infi, (6.17) 

u\an = p, 

in G{£1). To see this, we consider other representatives (a* J ), (6^), (pj, and (w e ) 
of [(a* J )], [(b\)], [(p e )], and [(w e )]. Then the proof of Proposition 13.81 clearly im- 
plies that 

N K 

- Di&Dfc) + Y. V ^Y> = Ve in n, (6.18) 

i,j = l i=l 

U~e\dU =Pe+ Va 

where r/ e E J\f(Q) and rj t is a net of functions satisfying (13.161) . But this implies 
that this choice of representatives also satisfies (16.171) in Q{Vt), so our solution 
[(u e )) is independent of the representatives used. 

□ 

This completes our proof of Theorem 14.11 We now conclude this article by giving a 
brief summary of everything that we have discussed. 
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7. Summary 

We began the paper with an example to motivate the Colombeau Algebra method 
of solving the semilinear problem (11.11) with potentially distributional data. In particu- 
lar, in Section [2] we proved the existence of a solution to an ill-posed critical exponent 
problem in Proposition 12.41 Our method of proving the existence of a solution to this 
problem consisted of mollifying the data of the original problem and solving a sequence 
of "approximate" problems with the smooth coefficients. We then obtained a sequence 
of solutions that yielded a convergent subsequence. The framework we used to obtain a 
solution to this problem was modeled on the more general Colombeau approach that we 
developed later in the paper, but required only basic elliptic PDE theory. Following this 
existence proof, we began to develop the more general Colombeau algebra framework. 
To this end, in Section 13.11 we introduced notation for Holder norms and stated two a 
priori estimates from [[51 that were made more precise by Mitrovic and Pilipovic in lfT2l . 
In Section 13.21 we then introduced the general framework for constructing Colombeau- 
type algebras and the Colombeau algebra Q{Vt) used in this paper. We then discussed a 
method used to embed the Schwartz distributions V'(VL) into Q{Vt). We used this embed- 
ding to analyze a problem of the form (11.11) with distributional coefficients. In particular, 
in Proposition |4.5| we determined explicit conditions under which we could solve a semi- 
linear problem of the form (11.11) with rough coefficients. Then we finished Section l3T2l by 
defining a class of semilinear operators on Q{Vl) in !3.4[ and we then defined the Dirichlet 
problem for these operators. 

Our main results for the Colombeau algebra framework were presented in Section 0] 
namely Theorem 14. 1[ which consists of an existence result for the semilinear prob- 
lem (14.21) in G(Vt). We then developed the necessary tools to analyze our semilinear 
problem in Section [51 We first determined a net of L°° bounds for positive solutions to 
our problem. Then, in Section [5T21 we showed that this net of L°° bounds is in fact a net 
of sub- and super-solutions contained in C, the ring of generalized constants described 
in Section [3T2l After developing our sub- and super-solutions, we proved Theorem 14.11 
in Section [61 We set up our problem in a manner similar to that used by Mitrovic and 
Pilipovic in [fT2l . However, our approach to solving our semilinear problem was distinct 
from theirs; we first determined a net of solutions (u e ) to the family of semilinear prob- 
lems (16.31) by using the method of sub-and super- solutions (Theorem 14.31) . and our net 
of sub- and super-solutions determined in Section 15.21 Once our net of solutions was 
determined, we then employed Theorems 13. II and our net of sub- and super-solutions to 
show that our net of solutions was contained in £m(Q)- 

In this article we have attempted to develop some basic tools to allow for a more 
general study of the Einstein constraint equations with distributional data. Our goal was 
to extend the current solution theory for scalar, critical exponent semilinear problems 
such as the Lichnerowicz equation, allowing for more irregular data than is currently 
covered by the existing solutions theories (cf. (71 [H for a summary of the known results 
for the CMC, near-CMC, and Far-CMC cases through 2009). As a next step, we hope to 
use the tools developed in this article to extend the near-CMC and Far-CMC existence 
framework for rough metrics developed in [[HI IH [10l [2]| to cover the rough data example 
studied by Maxwell in IfTTi . 
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